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ABSTRACT 

We present the first results from the largest deep extragalactic millimetre-wavelength 
survey undertaken to date. These results are derived from maps covering over 0.7 deg 2 , 
made at A = 1.1mm, using the AzTEC continuum camera mounted on the James 
Clerk Maxwell Telescope. The maps were made in the two fields originally targeted at 
A = 850 \im with the Submillimetre Common- User Bolometer Array (SCUBA) in the 
SHADES project, namely the Lockman Hole East (mapped to a depth of 0.9-1.3 mJy 
rms) and the Subaru XMM Deep Field (mapped to a depth of 1.0-1.7 mJy rms). The 
wealth of existing and forthcoming deep multi-frequency data in these two fields will 
allow the bright mm source population revealed by these new wide-area 1.1 mm images 
to be explored in detail in subsequent papers. Here we present the maps themselves, 
a catalogue of 114 high-significance sub- millimetre galaxy detections, and a thorough 
statistical analysis leading to the most robust determination to date of the 1.1mm 
source number counts. These new maps, covering an area ~ 3 times greater than the 
SCUBA SHADES maps, currently provide the largest sample of cosmological volumes 
of the high-redshift Universe in the mm or sub-mm. Through careful comparison, we 
find that both the COSMOS and GOODS North fields, also imaged with AzTEC, con- 
tain an excess of mm sources over the new 1.1mm source-count baseline established 
here. In particular, our new AzTEC/SHADES results indicate that very luminous 
high-redshift dust enshrouded starbursts (Si.imm > 3mJy) are 25-50 per cent less 
common than would have been inferred from these smaller surveys, thus highlight- 
ing the potential roles of cosmic variance and clustering in such measurements. We 
compare number count predictions from recent models of the evolving mm/sub- mm 
source population to these SMC surveys, which provide important constraints for the 
ongoing refinement of semi-analytic and hydrodynamical models of galaxy formation, 
and find that all available models over-predict the number of bright sub-millimetre 
galaxies found in this survey. 

Key words: surveys - galaxies: evolution - cosmology: miscellaneous - submillimetre 



1 INTRODUCTION 

In the last two decades, surveys in the far-infrared and sub- 
millimetre have revolutionised our understanding of galaxy 
evolution in the high-redshift Universe. These surveys, pri- 
maril y at waveleng ths ar ound 850-ixm (e.g. ISmail et al 



19971 : iHughes et all 1 19981 : iBarger et al.l 119981: IScott et al 
" l2003l ; iPope et all 120061; ICoppin et al 
. .Greve et all l2009h . 1100-nm (jLaurent et all 120051 : 
Scott et all |200S|; IPerera et all [20081), and 1200-n.m 



Greve et al.1 12004 iBertoldi etail 120071 : iGreve et all 120081 ') 

have shown that the contribution to the comoving IR 
energy density from sub-mm bright galaxies (SMGs) in- 
creases by approximately three orders of magnitude in go- 
ing from the local Universe to z <; 1, and that z J> 1 
SMGs are responsible for a significant portion the extra- 
galactic infrared ba ckground light (e.g. IScott et al. 1 120081 : 
ISerieant et al.ll2008l ). Initi al followup studies have identified 
optic al counterparts (e.g. iDve et all 2008; Clements et all 
2008 ) and measured redshifts (e.g. IChapman et al. 1 120051 : 



Aretxaga et al.1 120071 ) of many SMGs. Further studies have 



shown that SMGs harbour very high rates of star forma- 



Alexander et al. 20051: Kovacs et al.l 


20061: Coppin et all 


20081: iMenendez-Delmestre et al.1 12007 


: IPope et al.1 120081) 



and that at l east s ome are involved in ongoing mergers 
|Farrah et all 120021 : IChapman et~all 120031 ). implying that 
SMGs signpost massive galaxy assembly in the high red- 
shift Universe. Re view s of their prope r ties c an be found in 
iBlain et all (|2002h and lLonsdale et all l|2006l) . 



On a fundamental level, SMGs represent the efficient 
transformation of free baryons into stars and black holes. Re- 
cent work has therefore focused on understanding how SMGs 
relate to the cosmological evolution of the total and bary- 
onic mass density. Evidence suggests that this relationship 
is complex, with an intricate dependence on variables such 
as redshift, local environmental richness, and halo merger 
history. Accordingly, observations must find SMGs across a 
wide range of environments and redshifts. This has tradi- 
tionally proven difficult; SMGs are easy to find across wide 
redshift ranges due to the favourable k-correction at sub-mm 
wavelengths, but hard to find across wide ranges in environ- 
ment due to the inability of most sub-mm bolometer arrays 
to efficiently map large areas of sky. Constraints on SMG 
number counts have thus been limited by both sample size 
and cosmic variance, we have found relatively few of the 
brightest and rarest SMGs, and constraints on SMG clus- 
tering - an important tool in relating SMGs to the underly- 



ing dark matter distribution - 
|2004 Ivan Kampen et~aT1 120051 
also iFarrah et al.1 120061 



are weak 



Blain et al.1 
20091 . 



Chapman et al.1 !2UUiJI . see 
Magliocchetti et alj|2007l ). As a re- 
sult, many studies have adopted a two-pronged approach; 
using modest sized blank field surveys to constrain the prop- 
erties of the general SMG population, combined with tar- 
geted surveys of clusters to probe the properties of SMGs 
in the highest d e nsity regions (e.g 



Greve et al l 120071: iPriddev et al.1 120081 : lAustermann et al] 



Stevens et all 120031: 



20091 : iTamura et al.1 12009! ). Even this approach has draw- 



backs though, as it requires a pre-existing cluster catalogue 
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extending to redshifts significantly in excess of unity, where 
clusters are difficult to find. 

This situation has recently been imp roved by both the 
combined analysis of multiple surveys fe.g. lScott et al ]2006h 
and by the advent of larger area sub- mm surveys such a s 
the 850-u.m SCUBA/SHADES survey (|Coppin et al.ll2006T l. 
which mapped 0.2 deg 2 to depths of ogso ~ 2 mjy. In this pa- 
per, we present a further step forward in understanding the 
SMG population through the 1100-u.m AzTEC/SHADES 



survey, which covers 0.5 deg to depths of criio 



1 mjy 



and over 0.7 deg 2 in total. This survey dramatically im- 
proves our understanding of the 1100-Lim blank-field popu- 
latio n, with previous 1 100-u.m surveys being smaller in are a 
Ce.g. lPerera et al.ll200Sl) . shallower (e.g. lLaurent et al.ll2005l ). 
or c ontaining known biased regions in their surve y volumes 
(e.g. lAustermann et al.ll2009l ; iTamura et al]|2009h . 

This paper is organised as follows. Details of the AzTEC 
observations are given in Section[2] and we present the 1100- 
u.m maps and source catalogues in Section [3] Detailed con- 
straints on the 1100-u.m blank-field number counts are de- 
rived in Section [4] and we compare these results to those 
derived from other surveys and models in Section [S] Finally, 
our conclusions are summarised in Section [6] We assume 
a flat ACDM cosmology with Qm = 0.3, £Ia = 0.7, and 
Ho = 73 km s _1 Mpc -1 . This paper is the first in a series 
of papers using the AzTEC/SHADES maps and catalogues 
to study the sub-millimetre population of galaxies. 



2 OBSERVATIONS 

We have completed the SCUBA Half-D egree Extragalactic 
Survey (SHADES; iMortier et~aH I2005T ) by mapping over 
one-half square degree of sky using the AzTEC 1.1-mm 
camera (|Wilson et alj2008i ) mounted on the 15-metre James 
Clerk Maxwell Telescope (JCMT). The AzTEC/ JCMT sys- 
tem results in a Gaussian beam with FWHM « 18 arcsec. 
The SHADES survey is split between the Lockman Hole 
East field (LH; 10 h 52 m , +57°00') and the Subaru/XMM- 
Newton Deep Field (SXDF; 02 h 18 m , -05°00'). AzTEC 
has mapped over 0.7 deg 2 to 1.1 mm depths of 0.9-1.7 mjy 
between the LH and SXDF fields, including the central 
0.13 an d 0.11 deg 2 , respec tively, mapped by SCUBA at 
850 u.m (|Coppin et alj|2006h . All AzTEC/SHADES observa- 
tions were carried out between November 2005 and January 
2006, with over 180 hours of telescope time dedicated to this 
project, including all overheads. 

The AzTEC/SHADES observing strategy is similar 
to that used for other AzTEC blank-field surveys at the 
JCMT and is described in detail in pre vious publica- 
tions l|Scott et al.l l2008l ; [Perera et all l2008h . All observa- 
tions were made while scan ning the telescope in eleva- 
tion in a raster pattern ("see I Wilson et aU 12008) . Initially, 
the AzTEC/SHADES observations were made as small 
15 arcmin x 15 arcmin mosaic maps with scan speeds of 90- 
120 arcsec s _1 . Later observations were extended to cover an 
entire field in one continuous observation of size 35 arcmin x 
35 arcmin, which served to reduce observational overheads. 
Faster scan speeds of 180-220 arcsec s _1 were used for these 
longer scans, which increased the effective sensitivity of the 
observations due to the corresponding reduction in resid- 
ual atmospheric noise at the higher temporal frequencies 



|Wilson et alj|20 08f). After full reduction, the larger maps 
with faster scan-speeds have an observing efficiency of ~ 
150 per cent, relative to the smaller maps. In the end, 46 (63) 
mosaic maps and 65 (34) full maps were used to create the 
final LH (SXDF) map. These observations were performed 
over a wide range of atmospheric conditions and elevations. 
Figure [T] shows both the achieved mapping speeds and the 
cumulative distribution of observation time as a function of 
effective opacity at 225 GHz. 

Nightly overhead observations included focusing, load 
curves, beam maps and pointing observations , all of which 
are d escribed in the AzTEC instrument paper (jWilson et al.l 
2008). Pointing observations of bright point sources (typi- 
cally > 1 Jy) that lie near the science field being targeted 
were made every two hours. These measurements provide 
small corrections to the JCMT pointing model and are 
applied using a linear interpolation between the nearest 
pointing measurements taken before and after each science 
observation. Flux c alibration is performed as described in 
IWilson et al.l (120081 ) using the nightly load curves and beam 
maps of our primary calibration source, Uranus. The error in 
flux calibration is e stimated to be 6-13 per cent on an indi- 
vidual observation l|Wilson et al.l l2008). The actual error in 
the final co-added AzTEC/SHADES maps, which comprise 
observations spanning many nights and calibrations, will be 
smaller, assuming the calibration uncertainty is randomly 
distributed. These individual error estimates do not include 
the systematic 5 p er cent absolute unce rtainty in the flux 
density of Uranus (|Griffin fc Ortonlll993l ). 
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Figure 1. Achieved mapping speeds for the various 
AzTEC/SHADES observations. Mapping speeds do not in- 
clude overheads and are calculated as defined in IWilson et al.l 
(2008). The cumulative distributions of observation time go- 
ing into the final AzTEC/LH (92 hours) and AzTEC/SXDF 
(61 hours) maps are plotted as the solid and dashed curves, 
respectively, and normalised to an arbitrary value on the y-axis. 
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Figure 2. AzTEC Signal-to-Noise maps of the Lockman Hole East field and SXDF, each observed with an AzTEC/JCMT beamsize 
of FWHM PS 18arcsec. Maps are displayed at different scales. The most significant source candidates, as defined in the text and listed 
in Tables [T]{3] are circled. The observing pattern results in maps that are deepest in the centre, with noise increasing towards the 
edges. Reference contours show the distribution of la depth, in mjy. The maximal extent (trgso 6mjy) of the overlapping 850-um 
SCUBA/SHADES survey is depicted by the dotted curves, (a) AzTEC/LH S/N map with 0.45 deg 2 total displayed area. For analysis 
purposes, the survey area is trimmed to the region of relatively uniform coverage (thick 1.3 mjy contour), resulting in a 0.31 deg 2 map 
with noise le vels 0.9— 1.3 mjy. F or comparison in Section [5] the approximate size of the similar depth 1.1-mm AzTEC/GOODS-N survey 
(0.068 deg 2 ; iPerera et aljlioolj) is represented by the dash-dotted rectangle; however, these surveys do not actually overlap on the sky. 
(b) AzTEC/SXDF S/N map covering a total of 0.59 deg 2 , with the trimmed analysis region (thick 1.7mjy contour) covering 0.37deg 2 
with noise levels 1.0— 1.7 mjy. 



3 MAPS AND CATALOGUES 

In this section we describe the methods used to construct the 
1.1-mm maps and source catalogues. We test the astrometry 
and calibration of our maps against complementary radio 
data. We also describe expanded and improved methods for 
estimating and correcting for flux biases inherent to these 
surveys and test these estimates against simulations. 



3.1 Mapmaking 

The time streams of each observation are cleared of 
intermittent spikes (e.g. cosmic-ray events, instrumen- 
tal glitches) and have the dominant atmosph eric signals 
remov ed using the techniques described in IScott et al.l 
(2008). Each observation is then mapped to a 3arcsec 
x 3 arcsec grid in RA-Dec that is tangent to the ce- 
lestial sphere at (10 h 51 m 59 s , +57°21'43") for Lockman 
Hole and (02 h 18 m 01 s , -04°59'54") for SXDF. These are 
the same pixel sizes and tangent points used for the 
SCUBA/SHADES 850-um maps, allowing for straightfor- 
ward comparison of maps in upcoming SHADES publica- 
tions. All observations are then 'co-added' on the same grid 
to provide a weighted-average signal map and weight map 
for each field. 

In parallel, we pass a simulated point source (as de- 
fined through beam map observations) through the same 



algorithms to trace and record the effective PSF, or 'point 
source kernel', in our final maps. We also create five noise- 
only map realisations of each observation by jack-knifing 
(randomly multiplying by 1 or —1) each scan (5-15 sec- 
onds of data) of the time stream. This process works to 
remove any astronomical signal while preserving the domi- 
nant noise properties in the mapQ. The resulting jack-knifed 
maps are dominated by residual atmospheric contamination 
and detector noise. One-hundred fully co-added noise maps 
are then created by randomly selecting a noise realisation for 
each observation and calculating the weighted average in the 
same manner used to create the co-added signal maps. Be- 
cause the atmospheric contamination and the detector noise 
are uncorrelated amongst the full set of maps, the resulting 
co-added noise maps are, like the underlying noise in the 
signal map, extremely Gaussian. 

An optimal point source filter is applied to our maps 
utilising the information contained within the point source 
kernel and noise map realisations. The filtering techniques 
used are described in detail in previo us AzTEC publications 
ijScott et al.ll2008l : IPerera et al.ll2008T l. The resulting AzTEC 



1 Jack-knifing also removes confused astronomical signal, which 
can sometimes be considered a source of noise. However, confusion 
noise is not significant for maps of this depth and beamsize (see 
Section l4.2D , as confirmed through simulation. 
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Figure 3. Pixel S/N histograms of the LH and SXDF maps 
(thick histograms) and average of their respective noise map real- 
isations (thin histograms). The noise-only maps are well described 
as a Gaussian (smooth curve), while the signal maps are distorted 
by the presence of sources. Sources affect both the positive and 
negative flux distribution due to the zero-mean (AC-coupled) na- 
ture of these AzTEC maps. The distortions caused by sources are 
more apparent in LH due to the map's lower average noise level 
compared to SXDF (Figure 01. 



signal-to-noise maps of Lockman Hole and SXDF are shown 
in Fig. [2] The thick dashed contour of each map depicts the 
50 per cent coverage level, representing a uniformly-covered 
region that has a noise level within y/2 of that found in 
the deep central region of that map and beyond which the 
survey depth drops off sharply - a consequence of the par- 
ticular observation modes employed. The AzTEC/SHADES 
maps are trimmed at this 50 per cent coverage level - as de- 
fined by the co-added weight map described above - for all 
analysis in the following sections. The trimmed maps have 
total sizes of 0.31 deg 2 and 0.37 deg 2 for LH and SXDF, re- 
spectively, and correspond to depths of 0.9 < erih < 1.3 mjy 
and TO < owif < 1.7 mjy. The SXDF map is larger and 
shallower than that of Lockman Hole due to an observation 
script error that led to some individual maps being offset in 
declination. We continued to observe the resulting extended 
SXDF region after discovery of the error in order to max- 
imise the usefulness of our entire data set. 

The optimal filter is also applied to the co-added 
point source kernel and noise maps in order to provide the 
best model of the point source response and accurate esti- 
mates of the noise properties in our final map s. As shown 
for p r eviously publ i shed AzTEC/JCMT maps (|Scott et all 
l200Sl : [P crcra et all l2008h . these noise-only maps confirm 
a highly Gaussian nature of the underlying noise in the 
AzTEC/SHADES data (see Figure [3]). Accurate representa- 
tions of the filtered point source kernel and the noise prop- 
erties of the maps are critical components of the simulations 
described throughout this paper. 



3.2 Astrometry 

We have checked the astrometric accuracy of the 
AzTEC/SHADES maps by stacking (i.e. averaging) the 
AzTEC flux at the positions of r adio sources in th ese fields 
(techniques described in detail in lScott et al.ll2008l ). We use 
a re-reduction of the archiv al VLA 1.4 GHz c ontinuum data 
in the Lockman Hole field l)lbar et al.ll2009T l to generate a 
catalogue of radio sources in t he AzTEC/LH field and we 
utilise the 100 /iJy catalogue o f lSimpson etafl (|2006h in the 
AzTEC/SXDF field. These catalogues result in stacked de- 
tections of significance llcr and la for AzTEC/LH and 
AzTEC/SXDF, respectively. The stacked data are consis- 
tent with no systematic astrometric offset in either map 
with the possible exception of a small offset in declination, 
+2.9T1.3 arcsec, in the AzTEC/SXDF field. Due to the low 
significance and relatively small size of this potential offset, 
no correction is applied to the map. 

We can also constrain the random astrometric errors 
across the AzTEC maps by measuring the broadening of 
the stacke d signal compared to the AzTEC point-source re- 
sponse (see lScott et alfeOQct ). This broadening suggests that 
the random radial pointing error RMS is <r p <J 4 arcsec for 
both fields. Broadening of the stacked signal is also caused 
by pointing errors in individual AzTEC observations (final 
maps are co-additions of dozens of such observations) and 
clustering of the radio sources; therefore, <r p < 4 arcsec pro- 
vides a strong upper limit to the random astrometric errors 
in the AzTEC/SHADES maps. 



3.3 Calibration checks 

Flux density calibration was performed on a nightly ba- 
sis as described in Section [2] All steps in producing the 
AzTEC/SHADES maps were perform ed by the AzTEC 
instrument team l|Wilson et all [2008) and are identical 
to those employed for other publis hed AzTEC data sets 
( Scot t et al . 2008: lPerera et al.l |2008l). thus minimising sys- 
tematic differences between these AzTEC surveys. Flux cal- 
ibration is expected to be consistent across all regions of the 
final AzTEC/SHADES maps, with each point being sam- 
pled numerous times by each of a large set (^> 10) of ob- 
servations that span a wide range of atmospheric conditions 
and calibrations. 

We find no significant systematics between individ- 
ual observations, including between the small mosaic and 
large full-map observations. Noise properties are consistent 
across individual observations, with mapping speeds follow- 
ing the corre l ation with atmospheric opacity described in 
IWilson et al.l (l200Sf) . The metho ds used to remove atmo- 
spheric signal ( Perera et al.l[2008l ) produce consistent point 
source kernels (shape and amplitude) across all observations, 
with the resulting attenuation of point sources varying about 
the mean with an RMS of 3.6 per cent (attenuation corrected 
for in the final map based on the average reduction in flux 
density) . 

We also use the stacking analysis of Section 13.21 as a 
check of the relative calibrations between fields based on 
the average mm-wave flux of known radio populations. We 
find that the calibrations of the AzTEC/SHADES fields 
are consistent; the average 1.1-mm fluxes at the location 
of Si^ghz > 100 /^Jy radio sources are 0.59 ± 0.09 mjy and 



6 J.E. Austermann et al. 



0.50 ± 0.07 mjy for the LH and SXDF maps, respectively. 
If we remove radio sources that fall within 9 arcseconds of 
significant AzTEC 'sources' with S/N > 3.5 or S/N < -3.5 
from the stacking analysis, we find that the stacked aver- 
age 1.1-mm flux becomes 400 ± 100 ^Jy and 430 ± 70pJy 
for LH and SXDF, respectively. We also utilise the deeper 
radio catalogue available for the LH field to determine that 
the average 1.1-mm flux of millimetre-dim £i.4GHz > 66 /j.Jy 
radio sources is 528 ±71 fj,3y, which is consistent with similar 
stack s of the AzTEC/COSMOS (530 ± 87 ^Jy; IScott et all 
l200Sh and AzT EC/GOODS-N (439 ± 107 pJy; using the ra- 
dio catalogue of [Biggs fc Ivison|[2006l and the AzTEC map of 
iPerera et alfeoosh surveys. Assuming the various fields have 
similar radio populations, which may not be strictly true 
for fields with significan t structure (e.g. AzTEC/COSMOS; 
I Austermann et al]|2009l ). these tests show that the calibra- 
tion across various AzTEC surveys is consistent within the 
measurement errors of the stacks. Note that since the avail- 
able radio catalogue does not cover the entire AzTEC/LH 
field, we stack only on radio sources found in deep regions 
of the radio survey (<7i.4GHz < 13/iJy) to ensure a uniform 
sampling of Siaghz > 66 /iJy sources. 



3.4 Catalogues 

Candidate mm-wave sources are identified as local maxima 
in the optimally filtered maps that pass a chosen threshold. 
Although it is possible that some local maxima could be due 
to multiple neighbouring sources (relative to the FWHM 
w 18 arcsec beamsize) that blend into one peak, the low 
S/N of most detections prohibits deconvolution of poten- 
tial multi-source peaks. AzTEC maps are mean-subtracted 
(i.e. the background has a zero net contribution) and SMGs 
are expected to be sparse, with less than one source per 
AzTEC/JCMT beam down to < 0.1 mjy; therefore, source 
blending is expected to be rare for bright sources in the 
AzTEC/SHADES survey, unless the SMC population is sig- 
nificantly clustered on scales smaller than 18 arcseconds. 
Since very little is known about the SMC population on 
these scales, we caution the reader that the following analy- 
sis assumes no clustering at small scales. Interferometric ob- 
servations over the coming years with the SMA, CARMA, 
and later with ALMA will address the clustering question 
definitively. 

The most robust AzTEC/SHADES source candidates 
are given in Tables Q~r[3] Source candidates are listed in 
descending order of detected S/N. Centroid source posi- 
tions are determined using flux-squared weighting of the 
pixels within 9 arcsec (FWHM/2) of the local maxima. The 
AzTEC/SHADES survey detects 43 and 21 robust sources 
with S/N > 4 in the LH and SXDF fields, respectively. 
Additional significant detections, as defined in Section 13.61 
are also listed in the source tables and considered in the 
analysis of Section [4] Multi- wavelength analysis of sources 
found in the AzTEC/SHADES survey, including combined 
850 u.m/1100 u.m properties of sources within the overlap- 
ping AzTEC and SCUBA SHADES surveys (Negrello et al., 
in prep.), is deferred to future publications. 
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Figure 4. Bayesian posterior flux density (PFD; solid curve) 
compared to the intrinsic flux distribution recovered through 
simulation (histogram with la Poisson errors of simulation) for 
sources detected at the significance listed and at a noise level of 
<r m = 1.0 mjy. Also shown are the Bayesian predictions without 
the correction for the bias to peak locations in the map (dot- 
ted curves). The dashed vertical line represents the raw mea- 
sured flux, Sm, while the dashed curve represents the Gaussian 
probability distribution that might otherwise be assumed without 
flux boosting and/or false detection considerations. Negative flux 
probabilities are allowable through the zero-mean AzTEC point- 
source kernel and effectively represent the probability of a null 
detection. 

3.5 Flux corrections 

Sources discovered in these blind surveys experience two no- 
table flux biases, both of which cause the average measure- 
ment of the flux density of a detected source to be high rela- 
tive to its true intrinsic flux, Si. These flux biases can be very 
significant (~ 10-50 per cent for sources listed in Tables Q]- 
[3|, particularly for the low-significance detections that typify 
SMC surveys. Therefore, it is important to characterise and 
correct for these biases before fluxes and number counts can 
be compared to measurements at other wavelengths. Since 
these biases are a function of survey depth, these corrections 
are also necessary before detailed comparisons can be made 
with other 1.1mm surveys. 
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Table 1. AzTEC Lockman Hole (LH) candidates with S/N ^ 4. Additional robust AzTEC/LH sources with lower S/N values are 
listed in Table [2] The columns give: 1) AzTEC source name, including RA/DEC centroid position; 2) Nickname; 3) Signal-to-noise of 
the detection in the AzTEC map; 4) Measured 1100-um flux density and error; 5) Flux density and 68 per cent confidence interval, as 
defined in Section 13,61 after corrections for flux boosting and the bias to peak locations in the map; 6) Probability that the source will 
de-boost to S; < when assuming the AzTEC/SHADES Bayesian prior. 



Source 








Nickname 


S/N 


Shoo 
(measured) 
(mjy) 


Siioo 
(corrected) 
(mJy) 


ffSuoo < 0) 


AzTEC. 


.J105201. 


.98+574049. 


3 


AzLOCK. 


.1 


8. 


.2 


7. 


.4 + 0.9 




0.000 


AzTEC. 


J105206. 


.08+573622. 


6 


AzLOCK. 


.2 


8. 


,2 


7. 


.2 + 0.9 


6.4+o ; l 


0.000 


AzTEC. 


.J105257. 


.18+572105. 


9 


AzLOCK 


.3 


7. 


.5 


7. 


.2 + 1.0 


6-2+oi 


0.000 


AzTEC. 


.J105044. 


.47+573318. 


3 


AzLOCK. 


.4 


6. 


.7 


6. 


.2 + 0.9 


5.3+°j 


0.000 


AzTEC. 


.J105403. 


.76+572553. 


.7 


AzLOCK 


.5 


6. 


,4 


5. 


.9 + 0.9 


4.9+oi 


0.000 


AzTEC. 


.J105241. 


.89+573551. 


7 


AzLOCK. 


.6 


6 


,2 


5. 


.6 + 0.9 


4.8+i* 


0.000 


AzTEC. 


.J105203. 


.89+572700. 


5 


AzLOCK 


.7 


6. 


.0 


5. 


.7 + 0.9 


4.8+H 


0.000 


AzTEC. 


.J105201. 


.14+572443. 





AzLOCK. 


.8 


6. 


.0 


5. 


.6 + 0.9 


4.7+i : ° 


0.000 


AzTEC. 


.J105214. 


.22+573327. 


.4 


AzLOCK 


.9 


5. 


.7 


5. 


.0 + 0.9 


4.1 + o1 
— u.y 


0.000 


AzTEC. 


J105406. 


.44+573309. 


6 


AzLOCK. 


.10 


5. 


.5 


5. 


.1 + 0.9 




0.000 


AzTEC. 


J105130. 


.29+573807. 


2 


AzLOCK 


.11 


5. 


.4 


4. 


.8 + 0.9 


0+1.0 
9 


0.000 


AzTEC. 


J105217. 


.23+573501. 


.4 


AzLOCK. 


.12 


5. 


,4 


4. 


.7 + 0.9 


3.8+H 

— l.U 


0.000 


AzTEC. 


J105140. 


.64+574324. 


6 


AzLOCK 


.13 


5. 


.3 


5. 


.3 + 1.0 


4.1+1? 


0.000 


AzTEC. 


.J105220. 


.24+573955. 


1 


AzLOCK. 


,14 


5. 


.2 


4. 


.6 + 0.9 


— u.y 


0.000 


AzTEC. 


.J105256. 


.32+574227. 


5 


AzLOCK. 


.15 


5. 


.2 


i. 


.8 + 0.9 


3.7+tg 


0.000 


AzTEC. 


.J105341. 


.50+573215. 


9 


AzLOCK. 


.16 


5. 


,2 


4. 


.7 + 0.9 


3.7+i : ° 


0.000 


AzTEC. 


J105319. 


.47+572105. 


3 


AzLOCK 


.17 


5. 


.0 


4. 


.7 + 0.9 


3.6+i : ° 

— l.U 


0.001 


AzTEC. 


.J105225. 


.16+573836. 


7 


AzLOCK. 


.18 


4. 


.8 


4. 


.2 + 0.9 


3 2 +1 -° 


0.001 


AzTEC. 


J105129. 


.55+573649. 


2 


AzLOCK 


.19 


4. 


.8 


4. 


.3 + 0.9 


3 2 +1 '° 

— 0.9 


0.002 


AzTEC. 


.J105345. 


.53+571647. 





AzLOCK. 


.20 


4. 


,7 


4. 


.7 + 1.0 


3.4+J'l 


0.004 


AzTEC. 


J105131. 


.41+573134. 


1 


AzLOCK 


.21 


4. 


,7 


4. 


.1 + 0.9 


3-1+r'o 

— l.U 


0.003 


AzTEC. 


.J105256. 


.49+572356. 


7 


AzLOCK. 


.22 


4. 


,7 


4. 


.5 + 1.0 


3.2+H 

— l.U 


0.004 


AzTEC. 


.J105321. 


.96+571717. 


.8 


AzLOCK 


.23 


1 


.5 


4. 


.4 + 1.0 


3.1+i? 


0.007 


AzTEC. 


.J105238. 


.46+572436. 


.8 


AzLOCK. 


.24 


4. 


.5 


4. 


.3 + 0.9 


3.0±i : ? 


0.007 


AzTEC. 


J105107. 


.06+573442. 


2 


AzLOCK 


.25 


4 


4 


3. 


.9 ± 0.9 


2 7+ 1 ' 


0.008 


AzTEC. 


.J105059. 


.75+571636. 


7 


AzLOCK. 


.26 


4. 


,3 


4. 


.6 + 1.0 




0.016 


AzTEC. 


J105218. 


.64+571852. 


9 


AzLOCK 


.27 


4. 


3 


i. 


.3 + 1.0 


2.9+H 

— 1.1 


0.013 


AzTEC. 


.J105045. 


.11+573650. 


.4 


AzLOCK. 


.28 


4. 


3 


4. 


.1 ± 1.0 


2 7+ 1 ' 1 
z -'-l.l 


0.015 


AzTEC. 


J105123. 


.33+572200. 


.8 


AzLOCK 


.29 


4. 


.2 


i. 


.0 + 0.9 




0.016 


AzTEC. 


.J105238. 


.09+573003. 


.4 


AzLOCK. 


.30 


4. 


.2 


3. 


.8 + 0.9 




0.014 


AzTEC. 


.J105425. 


.31+573707. 


.8 


AzLOCK 


.31 


4. 


.2 


5. 


.2 + 1.2 


3 l +i:4 


0.038 


AzTEC. 


.J105041. 


.16+572129. 


6 


AzLOCK. 


.32 


4. 


,2 


4. 


.1 ± 1.0 


2 7+ 1 ' 1 
z -'-1.2 


0.020 


AzTEC. 


.J105245. 


.93+573121. 


2 


AzLOCK 


.33 


4. 


.2 


3. 


.8 + 0.9 


z -°-i.o 


0.016 


AzTEC. 


.J105238. 


.35+572324. 


.4 


AzLOCK. 


.34 


4. 


,1 


4. 


.0 + 1.0 


2.611;° 


0.023 


AzTEC. 


.J105355. 


.84+572954. 


.7 


AzLOCK. 


.35 


4. 


.1 


3. 


.7 + 0.9 


2.5+1;° 


0.021 


AzTEC. 


.J105349. 


.58+571604. 


3 


AzLOCK. 


.36 


4. 


.1 


4. 


.3 + 1.1 


9 7+1.2 
z -'-1.3 


0.032 


AzTEC. 


J105152. 


.72+571334. 


5 


AzLOCK 


.37 


4. 


.1 


4. 


.1 ± 1.0 


Z - D -1.3 


0.033 


AzTEC. 


J105116. 


.44+573209. 


9 


AzLOCK. 


.38 


4. 


.0 


3 


.5 + 0.9 


2 4+1-0 
z - 4 -l.l 


0.025 


AzTEC. 


.J105212. 


.26+571552. 


5 


AzLOCK 


.39 


1 


.0 


i. 


.0 + 1.0 


9 r+l.l 
z -°-1.3 


0.035 


AzTEC. 


.J105226. 


.58+573355. 





AzLOCK. 


,40 


4. 


.0 


3. 


.5 + 0.9 


9 Q + 1.0 

z -°-i.o 


0.026 


AzTEC. 


J105116. 


.34+574027. 


3 


AzLOCK 


.41 


1 


.0 


3. 


.7 + 0.9 


24 +i.i 

z '*-1.2 


0.036 


AzTEC. 


.J105058. 


.27+571842. 


.8 


AzLOCK. 


,42 


4. 





3 


.9 + 1.0 


24 +l.2 
z - 4 -1.2 


0.042 


AzTEC. 


.J105153. 


.10+572122. 


7 


AzLOCK. 


.43 


4. 


.0 


3. 


.8 + 1.0 


2 4 +l.i 
z -*-1.2 


0.039 



The primary flux bias in SMG surveys is commonly re- 
ferred to as 'flux boosting' and is due to the combination of 
a source density that increases sharply with decreasing flux 
and the blind nature of the su rvey (i.e. source s have previ- 
ously unknown positions); see iHogg fc Turner! |l998) for a 
full description of this effect. W e employ an advanced ver - 
sion of the Bayesian methods of ICoppin et alJ (|2005l . 120061 ) 



to correct for flux boosting and generate a full posterior 
flux density (PFD) probability distribution for each source 
candidate. The Bayesian approach requires a prior in the 
form of the assumed number density of sources projected 
on the sky (i.e. 'number coun ts') as a function o f flux. We 
use the iterative method of lAustermann et aD l)2009l ) to 
determine the most appropriate prior. We begin by using 
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Table 2. AzTEC/LH source candidates with S/N < 4. Columns as described in TablcQ] 



Shoo Shoo 
(measured) (corrected) 



Source 




Nickname 


S/N 


(mjy) 


(m 


Jy) 


^(Snoo < 


AzTEC 


J105241. 87+573406.1 


AzLOCK.44 


3 


9 


3.5 + 0.9 


2.3 


hl.O 
-1.1 


0.033 


AzTEC 


.J105154.82+573824.6 


AzLOCK.45 


3 


9 


3.5 + 0.9 


2.2 


-1.0 
-1.1 


0.035 


AzTEC 


.J105210.75+571433.8 


AzLOCK.46 


3 


9 


3.9 + 1.0 


2.4 


-1.2 
-1.3 


0.047 


AzTEC 


.J105306.80+573032.7 


AzLOCK.47 


3 


9 


3.6 + 0.9 


2.3 


-1.0 
-1.2 


0.037 


AzTEC 


J105431. 31+572543. 3 


AzLOCK.48 


3 


9 


4.0 + 1.0 


2.4 


-1.3 
-1.3 


0.052 


AzTEC 


.J105340.49+572755.0 


AzLOCK.49 


3 


9 


3.6 + 0.9 


2.2 


hl.l 
-1.1 


0.039 


AzTEC 


.J105205.59+572916.1 


AzLOCK.50 


3 


9 


3.5 + 0.9 


2.2 


-1.0 
-1.2 


0.042 


AzTEC 


.J105035.90+573332.1 


AzLOCK.51 


3 


9 


3.7 + 1.0 


2.2 


-1.2 
-1.2 


0.050 


AzTEC 


.J105206.79+574537.5 


AzLOCK.52 


3 


9 


4.2 + 1.1 


2.4 


-1.3 
-1.6 


0.070 


AzTEC 


J105435.20+572715.9 


AzLOCK.53 


3 


9 


4.0 + 1.0 


2.4 


-1.2 
-1.5 


0.062 


AzTEC 


.J105351. 57+572648.8 


AzLOCK.54 


3 


8 


3.5 + 0.9 


2.1 


hl.l 
-1.2 


0.050 


AzTEC 


.J105153.94+571034.3 


AzLOCK.55 


3 


8 


4.6 + 1.2 


2.4 


-1.3 
-2.0 


0.094 


AzTEC 


.J105203.84+572522.7 


AzLOCK.56 


3 


8 


3.6 + 0.9 


2.1 


hl.l 
-1.3 


0.055 


AzTEC 


.J105251. 38+572609.9 


AzLOCK.57 


3 


8 


3.6 + 0.9 


2.1 


-1.2 
-1.2 


0.056 


AzTEC 


.J105243.78+574042.6 


AzLOCK.58 


3 


8 


3.4 + 0.9 


2.1 


-1.1 
-1.2 


0.053 


AzTEC 


.J105044.92+573030.0 


AzLOCK.59 


3 


8 


3.4 + 0.9 


2.1 


-1.1 
-1.2 


0.054 


AzTEC 


.J105345.63+572645.8 


AzLOCK.60 


3 


8 


3.5 + 0.9 


2.1 


-1.1 
-1.2 


0.056 


AzTEC 


.J105257.19+572248.5 


AzLOCK.61 


3 


8 


3.7 + 1.0 


2.1 


-1.2 
-1.3 


0.063 


AzTEC 


.J105211. 61+573510. 7 


AzLOCK.62 


3 


8 


3.3 + 0.9 


2.0 


-1.0 
-1.2 


0.056 


AzTEC 


.J105406.14+572042.0 


AzLOCK.63 


3 


7 


3.7 + 1.0 


2.1 


-1.2 
-1.4 


0.074 


AzTEC 


.J105310.94+573435.6 


AzLOCK.64 


3 


7 


3.3 + 0.9 


2.0 


-1.1 
-1.2 


0.059 


AzTEC 


.J105258.39+573935.4 


AzLOCK.65 


3 


7 


3.4 + 0.9 


2.0 


-1.1 
-1.2 


0.061 


AzTEC 


.J105351.46+573058.2 


AzLOCK.66 


3 


7 


3.4 + 0.9 


2.0 


-1.1 
-1.2 


0.064 


AzTEC 


.J105045.33+572924.4 


AzLOCK.67 


3 


7 


3.4 + 0.9 


1.9 


-1.1 
-1.2 


0.065 


AzTEC 


.J105325.86+572247.3 


AzLOCK.68 


3 


7 


3.5 + 0.9 


2.0 


hl.l 
-1.3 


0.071 


AzTEC 


.J105059.74+573245.6 


AzLOCK.69 


3 


7 


3.3 + 0.9 


1.9 


-1.1 
-1.2 


0.064 


AzTEC 


J105121. 65+573333.6 


AzLOCK.70 


3 


7 


3.2 + 0.9 


1.9 


-1.0 
-1.2 


0.064 


AzTEC 


J105407.02+572957.7 


AzLOCK.71 


3 


7 


3.4 + 0.9 


1.9 


f 1.1 
-1.3 


0.071 


AzTEC 


.J105132.73+574022.1 


AzLOCK.72 


3 


7 


3.4 + 0.9 


1.9 


hl.l 
-1.2 


0.069 


AzTEC 


.J105157.08+574057.6 


AzLOCK.73 


3 


7 


3.3 + 0.9 


1.9 


fl.l 
-1.2 


0.068 


AzTEC 


.J105246.38+571742.5 


AzLOCK.74 


3 


7 


3.6 + 1.0 


1.9 


fl.2 
-1.4 


0.087 


AzTEC 


.J105309.72+571700.1 


AzLOCK.75 


3 


7 


3.6 + 1.0 


1.9 


fl.2 
-1.4 


0.087 


AzTEC 


.J105228.45+573258.0 


AzLOCK.76 


3 


7 


3.2 + 0.9 


1.9 


fl.O 
-1.2 


0.067 


AzTEC 


.J105148.13+574122.5 


AzLOCK.77 


3 


7 


3.3 + 0.9 


1.9 


fl.l 
-1.3 


0.073 


AzTEC 


.J105349.75+573352.4 


AzLOCK.78 


3 


7 


3.4 + 0.9 


1.9 


fl.l 
-1.3 


0.076 


AzTEC 


.J105232.60+571540.3 


AzLOCK.79 


3 


7 


3.6 + 1.0 


1.9 


fl.2 
-1.5 


0.088 


AzTEC 


.J105418.55+573447.5 


AzLOCK.80 


3 


7 


3.7 + 1.0 


1.9 


fl.l 
-1.5 


0.093 


AzTEC 


.J105321. 70+572308.3 


AzLOCK.81 


3 


6 


3.4 + 0.9 


1.9 


fl.l 
-1.4 


0.083 


AzTEC 


.J105136.91+573758.1 


AzLOCK.82 


3 


6 


3.2 + 0.9 


1.8 


fl.O 
-1.3 


0.079 


AzTEC 


.J105343.81+572543.6 


AzLOCK.83 


3 


6 


3.3 + 0.9 


1.8 


fl.O 
-1.4 


0.090 


AzTEC 


.J105230.53+572210.0 


AzLOCK.84 


3 


6 


3.4 + 1.0 


1.8 


fl.O 
-1.6 


0.099 


AzTEC 


.J105036.93+573228.9 


AzLOCK.85 


3 


6 


3.3 + 0.9 


1.8 


fl.O 
-1.5 


0.096 


AzTEC 


.J105037.18+572844.9 


AzLOCK.86 


3 


6 


3.3 + 0.9 


-i 7 +0.9 


0.099 



the SCUBA/SHADES (|Coppin et al.ll2006l ) 850-nm number 
counts, scaled to 1.1mm through an initial assumption of 
the 850/1100 urn flux ratio, as the initial prior. The prior is 
then iteratively adjusted using the empirical number counts 
of this survey (Section 14. 1[) , which quickly converges within 
a few iterations. As the widest-area deep millimetre survey 
to date, these iterative AzTEC/SHADES results provide the 
best 1.1-mm blank-field source number density prior avail- 
able. 

A second notable flux bias results from sources being 
defined as local maxima in the map. Since the position of 



the source is not independently known, nearby positive noise 
inevitably induces positional errors and this noise can com- 
bine with the off-centre beam-convolved flux of the source 
to outshine the true source being measured, thus resulting 
in an average positive flux bias in the local maximum that 
is taken as the measurement. The bias is independent of 
the aforementioned 'flux boosting' (the Bayesian prior is a 
noiseless calculation) and is instead a systematic of the ac- 
tual measurement, as opposed to an effect of the luminosity 
function bein g surveyed. This bi as to peaks (or 'noise gradi- 
ent bias', e.g. Ivison et al.ll2007l ) is minimised by optimally 
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Table 3. AzTEC/SXDF source candidates. Columns as described in Table [T] 



Shoo Shoo 
(measured) (corrected) 
Source Nickname S/N (mjy) (mjy) P(Snoo < 0) 



AzTEC.J021738. 


.52- 


-043330 


.3 


AzSXDF. 


1 


5. 


.2 


7. 


.4± 


1 


.4 


5 


.3 4 


-1.4 
-1 7 


0.002 


AzTEC.J021745. 


.76- 


-044747 


.8 


AzSXDF 


2 


4. 


.8 


5, 


,4± 


1 


.1 


4. 


.0+ 


-1.1 
-1 3 


0.003 


AzTEC.J021754. 


.97- 


-044723 


.9 


AzSXDF. 


3 


4. 


.8 


5, 


.3± 


1 


.1 


3. 


,8 4 


-1.2 

-1 2 


0.004 


AzTEC.J021831. 


.27- 


-043911 


.9 


AzSXDF 


.4 


4. 


.8 


6. 


,9± 


1 


.4 


4. 


.4-+ 


-1.7 
-1.6 


0.011 


AzTEC.J021742. 


.10- 


-045626 


.7 


AzSXDF. 


5 


4. 


.7 


5. 


.1 ± 


1 


.1 


3. 


6 4 


-1.2 
-1.2 


0.005 


AzTEC.J021842. 


.39- 


-045932 


.7 


AzSXDF 


6 


4. 


.6 


5. 


,8± 


1 


.3 


4. 


.0+ 


-1.3 
-1.6 


0.010 


AzTEC.J021655. 


.80- 


-044532 


.2 


AzSXDF. 


.7 


4. 


.6 


6. 


.3± 


1 


.4 


4. 


4 


-1.6 
-1.7 


0.019 


AzTEC.J021742. 


.13- 


-043135 


.6 


AzSXDF 


.8 


4. 


.5 


6. 


,4± 


1 


.4 


4. 


.0-+ 


-1.6 
-1.8 


0.025 


AzTEC.J021823. 


.10- 


-051136 


.7 


AzSXDF. 


9 


4. 


.3 


6. 


9± 


1 


.6 


3. 


,8 4 


-1.9 
-2.4 


0.061 


AzTEC.J021816. 


.07- 


-045512 


.2 


AzSXDF 


10 


4. 


.3 


4. 


,7± 


1 


.1 


3 


.1+ 


-1.2 
-1.3 


0.018 


AzTEC.J021708. 


.04- 


-045615 


.3 


AzSXDF. 


11 


4. 


.3 


5. 


.5± 


1 


.3 


3. 


,3 4 


-1.5 
-1.7 


0.039 


AzTEC.J021708. 


.03- 


-044256 


.8 


AzSXDF 


12 


4. 


.2 


5. 


,9± 


1 


.4 


3 


.3-+ 


-1.7 
-1.9 


0.053 


AzTEC.J021829. 


.13- 


-045448 


.2 


AzSXDF. 


13 


4. 


.2 


4. 


.4± 


1 


.1 


2. 


,8 4 


-1.2 
-1.4 


0.028 


AzTEC_J021740. 


.55- 
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Notes: a) Source is included in order to have a complete list of candidates with S/N 4, despite its relatively high null probability. 

where the fractional flux of the Gaussian beam (relative to 
maximum) is (5. Although this bias is relatively small in flux, 
it can have a strong effect on the Bayesian probability den- 
sities of low S/N detections (e.g. 50 per cent overestimate 
in probability that Si = 5 m for a S m = 3 ± 1 mjy mea- 
surement; see Fig- [4j . We note that the estimates provided 
by Equation [T] are significantly sma ller than the genera lised 
case provided by Equation B19 of llvison et al. | (20071) and 
are specifically tailored to AzTEC/ JCMT scanning observa- 
tions through simulation. Maps with a significantly different 
response to point sources (e.g. different beamsize or mapping 
strategy) may require a reevaluation of the functional form 
and parameter values of Equation [T] 

We correct for the bias to peak locations by subtract- 
ing fcpoak from the measured flux, 5* m , before calculating 
the Bayesian estimated PFD. The differences between the 
Bayesian PFDs with (solid curves) and without (dashed 
curves) this secondary bias correction can be seen in Fig. [4] 
for several S/N detection levels. Analysis of past surveys 
typically ignored this bias, but largely avoided its effects by 
restricting the analysis to only the most significant sources 
(e.g. S/N ^4). 

The de-boosted flux value listed for each source in Ta- 
bles[T}j3]is taken as the flux at the PFD local maximum near- 
est the detected flux, S m (see Fig. [4j| . Our improved estimate 



filtering the map for point sources, but can still be a signif- 
icant factor for low-significance sources. 

We characterise and quantify the bias to peak loca- 
tions through 10,000 simulations of the LH and SXDF 
maps. These simulated maps are generated by populating 
the noise-only maps with the flux-scaled point source kernel 
at random locations drawn from a uniform distribution and 
in accordance with a number counts distribution that is con- 
sistent with the final AzTEC/SHADES counts ( Section l4~Tj) . 
We generate simulated PFDs by cataloguing the input flux 
(Si) associated with each source measurement (5 m ,tr m ) re- 
covered in the simulated maps. These simulated PFDs are 
compared to the Bayesian estimate to characterise the re- 
maining bias (e.g. Fig. |3J), which comes primarily from the 
bias to peak locations. Through comparison of the PFDs 
over the flux range under investigation here (Si > 1 mjy) 
and for detections with S/N ^ 3 , we find that the aver- 
age flux bias incurred for an AzTEC/ JCMT measurement 
of (Sm,ffm) is well described by the equation 

&pcak(5m ' (Tm) = ^ eXP ( Z £f 1 ) (1) 

with a = 1 and f3 = 0.4. In this form, the bias is modelled 
as a effective independent noise elements that lie at a ra- 
dial distance from the true source that is equivalent to that 
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of the significant biases at work for low significance detec- 
tions leads to accurate PFDs down to at least S/N = 3, thus 
allowing us to utilise more of the maps' information when 
conducting the source-list driven number counts analysis of 
Section 0] 



3.6 Source robustness and false detections 

The effects of flux boosting (Section 13. 5 [1 make S/N a less 
than ideal measure of source robustness; therefore, we in- 
clude an estimate of the total probability that the source 
will de-boost to Si < 0, P(Suoo < 0), as determined from 
the integrated Bayesian PFD, in the source lists of Tables[T]- 
This provides a better metric than just S/N for the rel- 
ative robustness of the source detections, due to its depen- 
dence on both Sm and <T m , r ather than just the ratio S m /a m 
(see also lCoppin et al1l2006h . We have restricted Tables QJO 
to include only the most robust AzTEC/SHADES sources 
with P(Sixoo < 0) sC 0.1. The effective S/N, as a func- 
tion of <j m , of this 10 per cent 'null-threshold' is plotted in 
Fig. [5] We note that the absolute value of P(Snoo < 0) is 
highly sensitive to the Bayesian prior used. For example, 
if we instead assumed the result s of the relatively source- 
rich AzTEC/GOODS-N survey (jPerera et alJl200Sl ) as the 
Bayesian prior, the number of sources in AzTEC/LH pass- 
ing the 10 per cent null-threshold increases from 86 to 221. 
Therefore, it is important to consider the priors used when 
comparing the number of 'detections' in various surveys of 
this type. However, we note that the effect of the choice of 
prior on the resulting number counts (Section [4| is much 
less substantial, as the apparent change in the number of 
'detections' is largely counteracted by the survey complete- 
ness corresponding to the particular prior used. 

The number of false detections in a given source cat- 
alogue depends strongly on the chosen threshold for what 
is, and what is not, defined as a source. Due to the rela- 
tively large 18 arcsecond beamsize of AzTEC on the JCMT, 
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Figure 5. Effective S/N threshold as a function of noise 
level for the given null-threshold values when assuming the 
AzTEC/SHADES number counts as the Bayesian prior. These 
curves represent constant levels of robustness for an SMG detec- 
tion using AzTEC/ JCMT data. These values are unique to the 
particular Bayesian prior (number counts) assumed. 
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Figure 6. Ratio of the total number of detections to the 
number of significant noise-only peaks as a function of null- 
threshold for the 50 per cent coverage region of AzTEC /LH and 
AzTEC/SXDF. This provides an estimate of the relative number 
of false detections expected beyond a given source threshold, i.e. 
below a given null-threshold. 



the AzTEC/SHADES maps are expected to become signif- 
icantly 'full' of sources (on average one source per beam) 
when considering the expected high density of sources with 
1.1-mm fluxes < 0.1 mJy. Various estimates of the false 
detection ra t e of AzTEC/JCMT maps are explored in 
iPerera et al.l ((2008) , who conclude that the average num- 
ber of significant noise peaks in the jack-knifed noise-only 
maps provide a conservative overestimate of the number of 
false detections in the map (a consequence of true source sig- 
nal adding both positive and negative fluxes to the underly- 
ing noise distribution to our zero- mean maps). The ratio of 
number of sources in the signal map to the average number 
found in the corresponding noise-only maps is plotted as a 
function of null-threshold in Fig. [6] for the LH and SXDF 
fields. 



4 SMG NUMBER COUNTS 

In this section, we present the sky-projected densities of 
1.1-mm sources in the AzTEC/SHADES survey and the 
methods by which they are determined. These methods 
represent an extended and improved version of the algo- 
rith ms outlined in the S CUBA/SHADES number counts pa- 
per (|Coppin et al.|[200r3) . In Section r4.3l we provide paramet- 
ric fits to the number count results of the combined surveys. 
These number count results provide a useful measure of the 
SMG population through which we compare those found in 
other fields, at different wavelengths, and that predicted by 
various models and simulations (Section [5}. 
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4.1 Number counts: algorithm and results 

We calculate source number c ounts using the b o otstra p 
sampling methods outlined in lAustermann et all {2009), 
which are motivated by those used to determine the 
SCUBA/SHADES number counts jCoppin et alj|2006h . In 
this method, the catalogues of continuous sourc e PFDs are 
samp led at random and with replacement (e.g. IPress et al.l 
1992) in order to determine specific intrinsic fluxes for the 
sources in the catalogue. These samples are binned to pro- 
duce both differential and integral source counts as a func- 
tion of intrinsic flux. This sampling process is repeated 
100,000 times to provide sufficient sampling of the source 
count probability distribution. Sampling variance is injected 
by Poisson deviating the number of sources sampled in each 
of the 100,000 iterations around the actual number of sources 
detected in the map. Number count results are taken as the 
mean of each bin and the distribution across the iterations is 
used to characterise the associated uncertainty. The counts 
are then corrected for completeness, using estimates derived 
from simulation, and scaled for survey area. The resulting 
number counts are then taken as the new Bayesian prior 
and the entire process, including producing new catalogues 
of sources and their PFDs, is repeated in the iterative-prior 
process described in Section 13.51 For each iteration of the 
prior, the number counts are calculated for both the LH and 
SXDF surveys independently and also for the two surveys 
combined. The Bayesian prior chosen for the next iteration 
is always taken as the best fit to the combined result. This 
iterative procedure minimises our bias to the number counts 
assumed in the Bayesian calculations. 

Previous surveys us ing a si milar bootst ra pping 

technique llCoppin et all 120061 ; iPerera et all I2OO8I ; 
lAustermann et ah 120091 ) limited the source catalogue 
to those sources with negative flux probability of 
P(Si < 0) < 0.05, i.e. a null-threshold of 5 per cent. 
This null-threshold value was historically used to limit the 
number of false detections to a near negligible amount and 
to render the bias to peak locations relatively insignificant. 
However, the false detection probability is inherently 
accounted for in the bootstrap sampling method if accu- 
rate PFDs are used. As discussed in Section [3] our bias 
corrections result in PFDs that are accurate for all source 
candidates with S/N ^ 3, and possibly lower significance 
(currently untested below S/N = 3). The PFDs are partic- 
ularly accurate in the Si ^ 1 mjy flux range considered in 
this analysis. Since the traditional null-threshold of 5 per 
cent would limit the AzTEC/SHADES source candidate 
list to just those with S/N <; 4 (Fig. [5|, we explore fainter 
sources in the data set with the use of higher null-thresholds 
that incorporate a larger catalogue of source candidates in 
the derivation of source count densities. 

We derive combined AzTEC/SHADES number counts 
using null-thresholds of 5, 10, and 20 per cent. The 20 per 
cent threshold represents the lowest significance tested in 
our simulations (effective S/N <; 3) and safely avoids com- 
plications related to source confusion by keeping the density 
of detections sufficiently low. The results are consistent for 
all three threshold values tested and the variations between 
the results are, in general, much smaller than the formal 
68 per cent uncertainty of the number count estimates. We 
have verified through simulation that the use of the higher 



null-threshold values supplies additional data without intro- 
ducing any significant biases or systematics (Section 14. 2 [1 . 

The combined AzTEC/SHADES differential number 
counts using 5 per cent (open circles) and 20 per cent (filled 
circles) thresholds are compared in Fig. [7] The two results 
are nearly identical at high fluxes, but differ slightly in the 
lower flux bins. The variation at low flux is not surprising 
given that the 180 additional source candidates being con- 
sidered when using the softer 20 per cent threshold are all 
relatively low in flux, thus providing significantly more data 
in the lower flux bins. All AzTEC/SHADES number count 
uncertainties represent the 68 per cent confidence interval 
derived from the distribution of bootstrap iterations. All un- 
certainties assume a spatially random distribution of sources 
and, therefore, do not account for the effects of cosmic vari- 
ance/clustering. The differential number counts data points 
are strongly correlated, as described in Appendix |A"1 

Fig. [8] presents the integral source counts, N(> S), for 
both the LH and SXDF surveys using the 20 per cent null- 
threshold. Unlike the finite-bin differential counts measure- 
ments, the integral counts are threshold measurements (i.e. 
number of sources greater than flux, S) and can be de- 
rived at continuous values of flux. Therefore, the final com- 
bined AzTEC/SHADES results are depicted as a continuous 
68 per cent confidence region. The combined differential and 
integral number counts are also given in Table [4] with inte- 
gral counts listed at integer flux limits. 

4.2 Simulations and tests 

We test for biases and systematics in these techniques by 
applying the same number count algorithms to simulated 
maps of model source populations. Simulated maps are con- 
structed as described in Section 13.51 and we test a range 
of input model populations motivated by past 1.1-mm and 
850-u.m surveys. 

Our simulations show that the number counts estimates 
for any flux bin can be significantly biased towards the as- 
sumed value in the Bayesian prior, particularly if that bin is 
poorly sampled by the catalogue of source PFDs used to con- 
struct the number counts. We significantly reduce this bias 
in the lower flux bins by extending the sampled catalogue 

Table 4. AzTEC/SHADES differential and integral number 
counts, calculated as described in the text. The differential num- 
ber counts flux bins are l-mjy wide with effective bin centres 
(first column) weighted according to the assumed prior. Correla- 
tions amongst data points are described in Appendix IA1 
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Figure 7. Differential number counts for the combined 
AzTEC/SHADES survey in 1 mjy wide bins. The per cent value 
in parentheses represents the null-threshold used for each data 
set. The 5 per cent threshold results of AzTEC/SHADES have 
been artificially displaced to the left for clarity, but represent the 
same bins as the 20 per cent threshold results. Number counts 
for the AzTEC/COSMOS and AzTEC/GOODS-N surveys have 
been re-calculated using the final AzTEC/SHADES prior (solid 
cu rve), while applyin g the methods of this pa per to the data sets 
of IScott etaLl d20Qg| > and iPerera et al.l J2008T) . respectively, and 
are calculated for slightly different bins (i.e. not shifted) for im- 
proved clarity. Error bars represent 68 per cent confidence in- 
tervals. Bin centres are weighted by the assumed prior (solid 
curve). The thick horizontal dashed line represents the 'survey 
limit', defined as the source density that will Poisson deviate 
to zero sources 32.7 per cent of the time in a survey this size. 
The dot, dash, d ash-dot, and dash-dot-d o t-dot curves represen t 
the predictions oflRowan- Robinson! j2009l).lGranato et al.H2004) , 
Ivan Kampen et al. I (120051 1. and iBaugh et alj (120051) , respectively. 
The thick and thin dotted curves represent models with high- 
redshift formation limits of Zf = 4 and Zf = 5, respectively. 



to include fainter source candidates with P(Si < 0) values 
up to 20 per cent, thus providing more data in these oth- 
erwise poorly sampled bins. This is shown through example 
in Fig. [9] Although significant bias to the chosen prior can 
still be seen in the lowest flux bin (l-2mjy), this bin is still 
very sensitive to the 'true' population. Therefore, by itera- 
tively adjusting the prior based on the results (Section 13. 
we find that the bulk of this bias can be removed. This gen- 
eral result is also supported through full simulations with 
a precisely known input population. As expected, the re- 
sults based solely on the brightest source candidates (null- 
threshold of 5 per cent; open squares) are more severely 
biased by the assumed prior at low fluxes. 

The primary concerns when considering low-significance 
sources (e.g. S/N = 3.0) are: (a) false detections (noise 
peaks); and (b) source confusion (significant contribution 
from multiple sources in each measurement). However, false 
detections are inherently accounted for by having accurate 
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Figure 8. Integrated number counts for the AzTEC/LH and 
AzTEC/SXDF surveys with 68 per cent confidence intervals. 
Their combined constraint on the average blank-field number 
counts is shown as a continuous 68 per cent confidence region. 
A 3-parameter Schechter fit to the combined differential counts 
(Fig. is shown as the solid curve. Results of other 1.1— mm 
surveys are plotted for comparison, with the AzTEC/GOODS- 
N and AzTEC/COSMOS results being recalculated using the 
AzTEC/SHADES prior an d a 20 per cent nul l-threshold. The 
SCUBA/SHADES results dCoppin et alj l2006h are scaled to 
1.1 mm using the combined fits of Section l4.3l Discrete integrated 
number counts data points are calculated at varying flux values 
(i.e. not shifted) for increased clarity. Model predictions are plot- 
ted as described in Fig. [7] 

PFDs at the intrinsic fluxes being probed, and our simula- 
tions show that confusion does not play a significant role at 
fluxes Si > 1 mjy, based on an extrapolation of measured 
SMG number counts (e.g. ICoppin et al.l 120061 ; IPerera et al.l 
|200S| ; this paper) and the AzTEC/JCMT beamsize (FWHM 
= 18arcsec). Using the fitted results of Section \4. 31 the tra- 
ditional rule of thumb confus ion 'limit' o f 1 source per 30 
beams (f2 bea m = ^beam! e.g. iHogdbOQll ) is ~ 0.8 mjy for 
AzTEC/JCMT 1.1 mm data and is below the most likely in- 
trinsic fluxes of the individual sources considered here. Most 
importantly, our simulations find no significant systematics 
or biases between the input and output number counts of 
the constructed maps, thus confirming that neither of the 
above concerns present a problem for the AzTEC/SHADES 
results as given. 

4.3 Parametric fits 

For simulation and modelling of the SMG population, it 
is often useful to have a functional form for the number 
counts result. We fit the AzTEC/SHADES differential num- 
ber counts to the 3-parameter Schechter function 
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Figure 9. An example of the increase in sensitivity and reduction 
of bias in the low-flux number counts when using the 20 per cent 
null threshold in favour of the 5 per cent threshold. Here, the 
differential number counts of the AzTEC Lockman Hole survey 
have been re-calculated (symbols) when assuming a significantly 
different (and poor) prior that predicts a much lower number 
of faint sources (solid curve). For comparison, the dashed curve 
represents the prior used throughout the rest of this paper, which 
is very near to our best estimate of the true AzTEC/Lockman- 
Hole number counts. The various priors are displayed down to 
0.1 mjy to highlight the differences between the models and do 
not represent a measurement at these fluxes. The open symbols 
have been artificially displaced by +5 per cent along the x-axis 
for improved clarity. 



using Levenberg-Marquardt minimisation. We convert the 
normalisation parameter TV' to the more easily interpreted 
TV3mjy (the differential counts at Siioo = 3mJy) using the 
relation 



N: 



3mjy 



-3mJy/S' 



(3) 



The best-fit AzTEC/SHADES parameters are listed in 
Table The table also includes the results of a combined 
analysis of the currently available AzTEC 'blank-field' sur- 
veys AzTEC/SHADES and AzTEC/GOODS-N; however, 
the addition of the relatively small GOODS-N survey pro- 
vides only a slight increase in the constraint of the average 
SMG population. These results are relatively insensitive to 
the Schechter parameter a, which is strongly anti-correlated 
to, and somewhat degenerate with, the parameter S" in the 
flux range sampled (Si > lmjy). Therefore, we find that 
the AzTEC/SHADES number counts are nearly as well de- 
scribed by a Schechter function with the a parameter fixed 
to a reasonable v alue that is consiste nt with previous data 
sets (e.g. a = -2; ICoppin et al.ll2006t ). 

In previous incarnations of the bootstrap sampling 
method outlined in Section 14. 1| (ICoppin et al.l 120061 ; 
iPerera et all 120081 ; lAustermann et al.1 12009 ). formal fits to 
the differential number counts resulted in unrealistically low 
X 2 values, due to an underestimate of the correlations be- 
tween bins. We have now improved the algorithm for cal- 
culating the correlation matrix, which is described in Ap- 
pendix[5] However, the large correlations amongst the 1 mjy 



wide AzTEC/SHADES flux bins lead to a level of degener- 
acy that significantly complicates the application of typical 
fitting algorithms that incorporate the covariance matrix. 

We avoid such complications in the derivation of best-fit 
statistics by implementing a bootstrap sampling method of 
parameter uncertainty estimation that is similar to what is 
used in the error estimation for the individual number count 
data points (Section 14. 1|) . In this method, parameter space 
is explored by calculating best-fit parameters for each of the 
100,000 number count bootstrap iterations. Fig. [TO] shows 
the resulting parameter space for a two-parameter fit to the 
AzTEC/SHADES results using Equation [2] with Schechter 
parameter a fixed to a value of —2. Marginalised 68 per 
cent confidence intervals are used for the parameter uncer- 
tainties presented in Table [5] . We find that this alternative 
approach gives results that are comparable to that of formal 
fits, while providing a better characterisation of the true pa- 
rameter probability distributions by avoiding assumptions 
of Gaussian distributed uncertainty in the fitted parameter 
and number count errors. Since an explicit flux value is cho- 
sen for each source upon an individual iteration of the boot- 
strap (i.e. a single flux is chosen from the source's PFD), the 
number counts found by each realisation have flux bins that 
are effectively independent; therefore, this method provides 
a direct exploration of parameter space without necessitat- 
ing an explicit calculation of the bin-to-bin correlations that 
exist amongst the final averaged results of Table [4] 

Table [5] also includes the results of simultaneous fits 
to the AzTEC/SHA DES and published SCUBA/SHADES 
|Coppin et alJliooih results, where we have assumed the two 
surveys are sampling the same source population and that 
the number counts of the two bands are consistent within an 




Figure 10. Example distribution of best fit results of each it- 
eration of bootstrapped AzTEC/SHADES number counts. Fits 
are to Equation [2] with the Schechter parameter a fixed to —2. 
Contours represent the 68 and 95 per cent confidence regions. 
Vertical and horizontal lines represent the marginalised 68 per 
cent confidence intervals of S' and N^ m j y , respectively. 
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Table 5. Parametric fits to differential number counts of Table [4] using Equation [2] Uncertainties of AzTEC-only fits represent the 
marginalised 68 per cent confidence intervals derived from the distribution of bootstrap iterations (Section I4.3|l . The Schechter parameter 
a is given in Column 4, and is held constant for fit results given without a confidence interval. The quantity adust i s a f ree parameter 
representing the spectral index inferred by the simultaneous fit of the AzTEC/SHADES (1100-nm) and SCUBA/SHADES (850-iim) 
results, as described in the text. Uncertainties of the combined AzTEC and SCUBA fit are the formal la parameter errors when 
using Levenberg-Marquardt minimisation. The additional negative correction listed for a^ust represents an estimated correction for the 
systematic error induced by the SCUBA/SHADES choice of prior, as described in the text. All parameter values are for number counts 
at 1100-um. 
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average scaling of flux density. These fits are accomplished 
through the introduction of a free spectral index parameter, 
Qdust, which we have defined to reflect the average flux ratio 
between the two observing bands through the relation 



■S850 
Slioo 



As 



An 



(4) 



where Asso and Anoo represent the effective centre wave- 
lengths of AzTE C and SCUBA, respectively (see also 
iPerera et ai]|2008l ). The combined SHADES fit, assuming 
the nominal band centres of 850 u.m and 1100 |xm and using 
Levenberg-Marquardt minimisation, gives ad us t = 3.81 ± 
0.17 (5 , 85o/S'iioo = 2.67 ±0.12). The quoted uncertainties 
do not include systematic errors due to spectral differences 
between the SMGs and flux calibrators, which is expected 
to be smaller than the formal la uncertainty given, or any 
systematic calibration errors between the data sets (the la 
uncertainty of the formal fit is equivalent in size to a ~5 per 
cent systematic error in the measured flux ratios). For opti- 
cally thin thermal dust emission in the Ray leigh- Jeans limit 
(A 2> hc/kT), Qdust represents the dust emissivity index 
(adust = 2 + /3); however, the Rayleigh- Jeans approxima- 
tion is not strictly applicable at these wa velengths due to 



the expected temperatur e (Td ~ 35 K; e.g. Chapman et al.l 



120051 ; IKovacs et all 120061: ICoppin et al.ll2008l ) and redshift 
((z) ~ 2.2; e.g. Chapman et al. 2005]) of the typical SMC. 



We have tested for various other systematics between 
the two instruments' data sets by recalculating the AzTEC 
number counts under conditions and assumptions closely 
ma tching those existin g for the SCUBA/SHADES analy- 
sis l|Coppin et al.ll2006r i. Since the AzTEC/SHADES survey 
includes additional mapped area not covered by SCUBA, 
the comparison of the two surveys may be susceptible to 
cosmic variance on large scales (i.e. J> 0.1 deg 2 ). However, 
we find that there are no significant differences in the re- 
sults when restricting the AzTEC analysis to only those 
regions covered by SCUBA. We also find no significant dif- 
ferences when applying the same 5 per cent null-threshold 
(Section gT]) used in the SCUBA/SHADES analysis. The 
SCUBA analysis lacks a correction for the bias to peak map 



locations (Section 13. 5[l . however, their use of the conserva- 
tive 5 per cent null-threshold should keep this bias relatively 
small and it is expected to have no significant effect on the 
resultant number counts. 

Finally, we note that the SCUBA analysis uses an ex- 
ternal Bayesian prior that was based on 850-u.m r esults of 
the Hubble Deep Field North (|Borvs et al.ll2003l ). as op- 
posed to the self-consistent iterative prior used in this pa- 
per. This prior represents a slight overabundance of bright 
SMGs when compared to the SCUBA/SHADES number 
counts, probably resulting in a small, but systematic, over- 
estimate of the SCUBA/SHADES counts. Although we can- 
not use the exact same prior as SCUBA without inher- 
ently assuming an 850 (i.m/1100 itm scaling relation (e.g. a 
value of ctdust), we can adopt a similar prior that assumes 
the results of an 1100-iim survey of the sam e approximate 
field (AzTEC /GOODS-N; [Perera et all [20081 ). We recalcu- 
late the AzTEC/SHADES counts with this prior, as well 
as other matched systematics (5 per cent null-threshold, no 
correction for peak bias) to re-estimate adust- These val- 
ues are compared to the previous fits to determine the sys- 
tematic error estimates given in Table [5] and result in a 
final corrected scaling index of cVd us t ~ 3.6 ± 0.2. As dis- 
cussed in Section [S] this value of adust is significantly larger 
than that inferred b y current measuremen ts of the SMG 
redshift distribution (Chapman et al. 2005) and the SEDs 
of local starbursts (|Dunne fc Ealesl 120011). as well as di- 
rect measurements of SMG flux ratio s (IKovacs et aL 2006; 
ICoppin et all l2008i ; iGreve et al.ll2008l : IChapin et all 120091 ). 
The lower Ssso /Shoo flux ratios detected through a direct 
comparison of SMGs det ected in the GOOD S-N field by 
both AzTEC and SCUBA (|Chapin et al.ll2009l ) indicate that 
our relatively high inferred flux ratio may be limited to the 
comparison of number counts and not due to systematic cal- 
ibration errors between the two instruments. This suggests 
that the differences between the AzTEC and SCUBA num- 
ber counts analyses and/or our assumptions of the source 
population (i.e. uniform flux ratio and the two wavebands 
track the same SMG population) lead to significant system- 
atic errors in the inferred flux ratio at this level of sensitivity. 
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Analysis of the 850 urn/ 11 00 um flux density ratios of indi- 
vidual SHADES sources is deferred to Negrello et al. (in 
prep.). 



5 DISCUSSION 

5.1 Comparison of 1.1mm surveys 

To appreciate the contribution of AzTEC/SHADES to our 
understanding of the SMG population, it must be compared 
to previous SMG surveys. The AzTEC/SHADES integral 
number counts are in strong disagreement with the para- 
metric results derived from a fluctuation analysis of the 
1.1-m m BOLOCAM Lockm an Hole survey (dashed line in 
Fig. 151 iMalonev et al.|[2005l ). The BOLOCAM survey is sig- 
nificantly smaller and shallower than this AzTEC/SHADES 
survey and consequently contains fewer sources and is more 
susceptible to sample variance. Furthermore, the BOLO- 
CAM fluctuation analysis is likely to be skewed by their 
requirement that the source population be well described 
by a single power law, which diverges at zero flux and has 
since been shown to poorly describe the SMG population 
over a wide range of flux densities (e.g. IScott et alj 120061 ; 
ICoppin et all 120061 ; this data set). Therefore, the BOLO- 
CAM/LH single power-law result may represent a compro- 
mise between the relatively steep drop in SMG number 
counts at high flux density and the inevitably more mod- 
erate slope at the faint end. 

Fig. [8j also shows the integral number counts for the 
individual AzTEC/LH (filled circles) and AzTEC/SXDF 
(open circles) fields at specific flux density limits. The two 
fields' number counts are consistent within their respective 
uncertainties; however, the overall trend suggests that the 
AzTEC/LH field is rich in bright (Shoo ii 4mJy) sources 
relative to AzTEC/SXDF. This difference of bright AzTEC 
source counts is consistent with the differences seen between 
the regions of LH an d SXDF surveyed by SCUBA at 850 um 
jCoppin et al.ll2006l ). 

The effects of cosmic variance appear more prevalent 
when comparing the results of this survey to the 0.15 deg 2 
AzTEC/COSMOS survey (|Scott et all 120081 ) . which tar- 
geted a region with significant structur e, as traced by the op - 
tical/IR galaxy population at z <J 1.1 (|Scoville et al.ll2007f) . 
The average blank-field number counts of AzTEC/SHADES 
confirm the significant over-density of bright 1.1-mm 
so urces in the AzT E C/CO SMOS region first reported 
in lAustermann et ail ((2009) , who conclude that the ob- 
served over-density is probably due to gravitational lens- 
ing by foreground (z <, 1.1) structure (comparisons be- 
tween SHADES sources and other populations/structure 
will be explored in future SHADES papers). We have re- 
calculated the AzTEC/COSMOS number counts using the 
AzTEC/SHADES prior and a 20 per cent null-threshold 
(Fig. |8j, affirming that the AzTEC/COSMOS over-density 
is significant regardless of the chosen prior. 

We similarly find that the AzTEC/GOODS-N re- 
gion is relatively rich in 1.1-mm sources compared to 
the much larger AzTEC/SHADES survey. This is consis- 
tent with the relative abundances found in the com para- 
ble 850-u. m surveys of GOOD S-N l|Borvs et al.l I2003T ) and 
SHADES (jCoppin et al.ll2006l ). The higher number counts 



of AzTEC/GOODS-N may be due to sample and/or cosmic 
variance on the scale of the GOODS-N map (0.068 deg 2 to 
~ 1 mjy) , which can be exemplified by moving a box the size 
of AzTEC/GOODS-N to different locations within the well- 
covered, and similar depth, regions of the AzTEC/LH map 
(e.g. dash-dotted rectangle in Fig. 2(a) I. This simple exercise 



shows that the total number of source candidates within the 
GOODS-N sized box can change by a factor of ~ 2 for any 
of the source definitions explored here (i.e. null-threshold 
5-20 per cent). The relatively high number of bright SMGs 
found in GOODS-N may be due, in part, t o potential high- 
redshift structures in th e GOODS-N field (|Chapman et al.l 
l2009l ; lDaddi et alj|2009h . 

To better quantify the e mpirical variations across fields, 
we turn to the formalism of Ef stathiou et al.l fl990). Tak- 
ing the three 'blank-fields' considered here (AzTEC/LH, 
AzTEC/SXDF, AzTEC/GOO DS-N) as independent 'c ells' 
and using Equations 9 & 5 from lEfstathiou et al.1 l|l990T ). we 
can calculate a 2 and Var(<7 2 ), where a 2 is defined by Equa- 
tion lb of the same paper and represents the variation from 
cell to cell, or in this case field to field. The variations in 
field size and completeness are taken into account as per the 
'counts in cells' formalism of Efstathiou (1990). For Var(<r 2 ) 
we assume zero variance {a 2 — 0) such that the significance 
of the measurement is essentially the 'detection' significance 
of some field to field variation. 

Number count data at the lowest fluxes, Shoo < 2 mjy, 
are omitted from this analysis due to their sensitivity to the 
assumed Bayesian prior, which is held constant for all fields 
(Section [3]). Any bias to the prior would result in measured 
field-to-field variations that are systematically lower than 
that expected from Poisson statistics. This bias is relatively 
small for the remaining data (Shoo > 2mJy) and, to the ex- 
tent that it is present, would act to make the fields' measured 
number counts systematically less varied and our empirical 
variance measurements conservatively low. We combine the 
counts such that we test the variations in two relatively well- 
sampled flux bins: 2-4 mjy, and > 4 mjy. 

Table [6] contains the results of this variance analysis. 
For the AzTEC 'blank-fields' alone, no significant detection 
of variance is found. The 2-4 mjy flux bin has a a 2 that is 
strongly negative, indicating that the measured variance of 
that bin is significantly less than that expected from Pois- 
son statistics, while the > 4mJy flux bin has a variance 
consistent with a 2 — 0. The former may be a consequence 
of the fact that AzTEC/LH and AzTEC/SXDF have coin- 
cidentally similar number counts compared to their formal 
uncertainty in the 2-4 mjy range. 

The results when including the AzTEC/COSMOS data 
as an additional cell in the analysis are also given in Table 
[6] It can be seen that for the brightest sources Shoo > 4 
mjy some variance is detected at the 3.8<r level. This fur- 
ther confirms the significant over-density of bright 1.1-mm 
sources in the AzTEC/COSMOS region. 

Interestingly, with the exception of the > 4 mjy sources 
across all 4 fields, these results are fairly consistent with 
what would be expected from consideration of the expected 
form of the correlation function. The expected variance can 
be calculated by integration of the correlation function 



j f t{r)dVdV, 
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Table 6. Statistical significance of the variations in the num- 
ber counts between the three fields considered here (LH, SXDF, 
GOODS-N), as well as the COSMOS field 



Flux a 2 A(cr 2 ) 
(mjy) 


Significance 

(A) 


N' 

(deg~ 2 ) 


LH, SXDF, GOODS-N only 


2-4 -8.9 x 1CT 3 2.3X1CT 2 
>4 7.9xl0" 2 8.2xl0" 2 


0.96 


464.5 
50.3 


LH, SXDF, GOODS-N and COSMOS 


2-4 l.lxlCT 2 2.3X1CT 2 
>4 2.4xl0 _1 6.4xl0~ 2 


0.48 
3.8 


489.9 
62.2 



Notes: Negative a 2 indicates the measured variance is less than 
expected from Poisson statistics. 

where the integral is calculated over a volume, V, defined 
as a truncated cone of solid angle SI over the redshift range 
2 < z < 3. We assume a correlation function of the form 
f(r) = (r/ro)~' r , with t~o = 5h _1 Mpc and 7 = 1.8, which are 
typical for local galaxy populations. Assuming field sizes in 
the range 0.07-0.37 deg 2 . gives a predicted variance of be- 
tween 0.017 and 0.008, respectively. The measured variance 
for the number density of Shoo > 4mJy sources across all 
4 fields is in significant excess of that predicted under the 
above assumptions. 

However it is worth noting that the quoted errors on 
the measured a 2 assume no clustering, and are therefore 
underestimates of the true measurement error. In addition 
it is known that the COSMOS field contains a significan t 
over-density of 1.1-mm sources (| Austermann et all [2009). 
Taken together it is clear that the COSMOS field is simply 
a highly unusual example, and the volumes probed by these 
surveys are not great enough to detect a clustering signal 
of bright SMGs through the comparison of number counts 
alone. 



nent being ~ 20 K. For a population of SMGs residing at the 
typical redshift of z ~ 2, the observed 850-(i.m SCUBA band 
is sampling a rest-frame wavelength of ~ 280 \im. To produce 
a similar Qdust to the SLUGS galaxies in the local Universe, a 
much hotter temperature is required for the SED (T ^ 50 K 
with f3 — 2). Alternatively, these SMGs could have simi- 
lar SEDs to the local galaxies but reside at lower redshifts 
(z <J 1), although this would be inconsistent with measured 
SMG redshift distributions (|Chapman et al.l 120051 ). The in- 
ferred sub-m m/mm flux ratio i s high compared to the model 
predictions of lSwinbank et all (|2008l ) an d at odds with mea- 
surements of the fl ux ratios ffasn /Sum jKovacs et aljfeood : 
ICoppin et all 120081), SWSii oo jChapin et al.l l2009h . and 
S850/S1200 ( Greve et alj|2008l ). which are all more consis- 
tent with the SLUGS SEDs for z ~ 2. In addition, the ex- 
iste nce of a populatio n of sub-millimetre drop-outs (SDOs; 
e.g. iGreve et al.ll2008h - sources with a combination of high- 
redshift and/or low dust temperature such that the 850-(i.m 
band samples near, or shortward, of the peak emission - 
would act to lower the average value adust for millimetre 
detected sources. 

It thus appears that our estimate of ad ust is systemati- 
cally large given the expectation that j3 2 |Dunne fe Ealesl 
l200ll . and references therein) and that the majority of our 
sources are unlikely to be fully in the Rayleigh- Jeans limit. 
This bias may be indicative of further systematics in the 
SCUBA/SHADES choice of prior (Section [O]), potentially 
insufficient de-boosting of low S/N SCUBA detections (as 
suggested in a direct comparison of sources det ected by both 
AzTE C and SCUBA in the GOODS-N field; IChapin et ail 
2009), or that selection bias somehow results in a system- 
atic increase in the value of adust inferred from SMG num- 
ber counts when assuming 850 \im and 1100 U-m sample the 
same approximate source population. A straight comparison 
of the AzTEC and SCUBA SHADES maps (Negrello et al., 
in prep.) will provide a more direct measure of ad us t that 
is based on individual sources and fluxes in the maps and 
search for evidence of SDOs in the SHADES fields. 



5.2 Comparison to 850— urn counts 

As shown in Fig. @ the SCUBA 850-^m and AzTEC 1100- 
11m SHADES counts are consistent within a uniform scal- 
ing of flux density (Section 14. 3p . Under the assumption that 
AzTEC and SCUBA are sampling the same source popu- 
lation (ignoring selection effects), adust represents a power 
law approximation to the average redshifted SMG SED 
at observed wavelengths of ~ 1 mm. The relatively steep 
850 urn/1100 u.m spectral index derived from the SMG pop- 
ulations of SHADES (after an approximate correction of sys- 
tematics due to chosen priors; see Section|4]3]), adust ~ 3.6± 
0.2, is roughly consistent with the 450 urn/850 [im spectral 
index, adust ~ 3.6-3.7, found by the SCUBA Local Universe 
Galaxy S urvey (SLUGS) of IR bright galaxies in the local 
Universe l|Dunne fc Ealesl l200ll ) after correcting for an av- 
erage CO (3-2) contamination of 25 per cent in the 850-(i.m 
band at z = (Se aquist et al.l I2004T ). However, the 450 \im 
measurements from SLUGS are already shortward of the 
Rayleigh- Jeans limit (but longward of the peak) in the local 
Universe; particularly as local galaxy SEDs require two or 
more dust temperature components, with the cooler compo- 



5.3 Predictions from models 

Finally, we compare the AzTEC/SHADES number counts 
to those predicted at 1100 um by various IR/sub-mm forma- 
tion and evolution models in figs. [7] and [8] The predict i ons o f 
the IR/sub-mm evolution models of iRowan- Robinson! (|2009l ) 
are shown for high-redshift formation limits of Zf = 4 and 
z f = 5. The AzTEC/SHADES number counts agree with 
the Z{ = 4 model at fluxes Snoo ^> 4mJy, but are system- 
atically lower than the predictions at higher fluxes. A semi- 
analytical model for the joint formation and evolution o f 
spheroids and QSOs (|Granato et al.ll2004l ; Isilva et al.ll2005l ) 
predicts very similar number counts at 1100 urn. These mod- 
els are in better overall agreement with the high flux number 
counts seen in the AzTEC/COSMOS and AzTEC/GOODS- 
N fields; however, those fields have significant biases and/or 
limitations, as discussed above. Also compared are the 
counts pr edicted by th e sem i-analytical galaxy formation 
model of iBaugh et al] (2005, see also lLacev et al. I 120081 ; 
ISwinbank et alj " 20081 ). which systematically over-predicts 
the number of sources seen in AzTEC/SHADES by a factor 
of 3-4 at all measured fluxes. Finally, we compare our results 
to the early predictions for SHADES (|van Kampen et al.l 
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120051 ) - models constr ained to the SC UBA 8-mJy (i.e. 
Shoo ^ 3mJy) survey (|Scott et alj|2002h - which forecast 
a shallower slope in the number counts than seen in the 
AzTEC/SHADES fields. Assuming the bright sources are 
uniformly distributed across the sky, the AzTEC/SHADES 
survey suggests that all of these models significantly over- 
predict the number of intrinsically bright SMGs. If, instead, 
these relatively rare sources are strongly clustered, the true 
all-sky average number density of the brightest SMGs could 
be higher (or lower) than indicated by this survey, poten- 
tially bridging the gap between model and observation. 



6 CONCLUSIONS 

AzTEC/SHADES is the largest extragalactic millimetre- 
wave survey to date, with over 0.7 deg 2 mapped to depths 
of 0.9 < oi.i < 1.7mJy. This survey, split between 
the Subaru/XMM-Newton Deep Field and the Lockman 
Hole, provides over 100 significant individual detections at 
1.1mm, with most representing newly discovered mm- wave 
sources. These maps also provide information on the fainter 
SMG population through the signature of numerous dimmer 
sources that are partially buried in the noise. 

Combined with our improved methods for number 
count estimates, AzTEC/SHADES provides the tightest 
available constraints on the average SMG population in 
the flux range lmjy <j Shoo lOmJy. In particular, 
the AzTEC/SHADES results represent a significant ad- 
vance in our knowledge of the blank-field population at 
1.1mm, showing that there are significantly lower densi- 
ties of bright SMGs than that suggested by smaller 1.1- 
mm surveys published previously. An accurate understand- 
ing of the average SMG population is critical for compar- 
isons to source counts found in biased and/or over-dense 
regions. The AzTEC/SHADES blank-field counts confirm 
the over-density of Sh op > 2 mjy sources fou nd in the 
AzTEC/COSMOS field jAustermann et alj|2009l ) and show 
that the GOODS-N field is also relatively rich in bright 
SMGs, thus suggesting that cosmic variance can signif- 
icantly affect the observed number density of SMGs in 
mass-biased regions (AzTEC/COSMOS) and/or on rela- 
tively small scales (AzTEC/GOODS-N; 0.068 deg 2 ). We find 
that the variance in number counts seen across the four avail- 
able AzTEC/JCMT survey fields (LH, SXDF, COSMOS, 
GOODS-N) is significantly larger than that expected from 
Poisson statistics alone, particularly at Shoo > 4mJy, thus 
suggesting that bright SMGs may be strongly clustered. 

The AzTEC/SHADES results are consistent with the 
predictions of the formation and evolution models of 
iGranato et al. I (I2004T ) and the zi ~ 4 evolution models 
of iRowan- Robinson! (|2009T ) for blank-field 1.1-mm source 
counts at Shoo Si 4mJy; however, these models systemati- 
cally over-predict the number of AzTEC/SHADES sources 
seen at higher fluxes, although the relative scarcity and po- 
tential clustering of bright sources leaves even this unprece- 
dentedly large SMG survey susceptible to the effects of cos- 
mic variance. A truly unambiguous characterisation of the 
Shoo il 6mJy SMG population will require significantly 
larger-area surveys at (sub-)mm wavelengths, such as those 
expected to be conducted in the coming year(s) by SCUBA- 



2 on the JCMT and AzTEC when mounted on the Large 
Millimeter Telescope (LMT). 

We find that the SCUBA/SHADES and 
AzTEC/SHADES number counts are consistent within a 
uniform scaling of flux density. Assuming that the 850 [im 
and 1100 u.m wavebands sample the same underlying source 
population, this scaling corresponds to an average source 
flux ratio of Ssso/Snoo ~ 2.5 ± 0.1, once corrected for 
known systematics between the data sets. This ratio is 
significantly larger than that expected for the high-redshift 
SMG population and we find that the systematics induced 
by small differences in the number count analyses of the 
two surveys and the assumption of a uniformly scalable 
flux density limit the robustness of the inferred flux ratio. 
The S850/S1100 flux ratio is explored further in a direct 
comparison of individual sources lying in the overlapping 
regions of the SCUBA and AzTEC surveys (Negrello et al., 
in prep.). 
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APPENDIX A: CORRELATION MATRIX 

The bootstrap sampling method of Section 14.11 induces sig- 
nificant correlation between the final averaged differential 
number count bins (e.g. Table [3} through discrete sampling 
(and consequent binning) of continuous PFDs that have sig- 
nificant probability on scales comparable to, or larger than, 
the bin size (lmjy in this p aper). Previous inca rnations of 
this sampling method (e.g. ICoppin et alj |2006l ) estimated 
covariance and correlation matrices directly from the vari- 
ation in number count results seen across the iterations of 
the bootstrap. This sampling method collapses each source's 
probability distribution (PFD; e.g. Fig. 3J) to a single flux 
upon each iteration, which acts to hide significant correla- 
tion amongst the final binned results by throwing away much 
of the cross-bin information contained within the PFD. This 
resulted in severely underestimated bin-to-bin correlations 
amongst the differential number counts data, as was evi- 
denced by the unrea l istically low v 2 val u es of formal fitting 
dCoppin et al.l 120061 ; iPerera et aTll2008l ; lAustermann et al.l 
2009). 

We now present an alternative method of calculating 
the correlation matrix which better captures these correla- 
tions amongst the final differential number count bins. We 
begin by integrating the PFD of each source over the span 
of each flux bin. These binned probabilities can be summed 
over all sources to provide a number counts estimate that 
matches the final averaged results of the full bootstrapping 
method, but without the robust uncertainty estimates that 
the bootstrap method is designed to provide. We apply this 
alternative number counts estimate to each of the 100,000 
unique catalogues produced by the bootstrap. The Poisson 
deviation and replacement sampling used to produce each 
catalogue (Section 2TJ act to perturb this new estimate of 
the number counts around the most likely values. This col- 
lection of perturbed number counts is then used to estimate 
the correlation between the differential number count bins. 

We present the resulting correlation matrix for the 
AzTEC/SHADES differential results (Section |4j in Ta- 
ble \Kl\ These correlations apply directly to the differential 
results provided in Tableland Fig. [7] (20 per cent thresh- 
old counts). The last row of Table |A"T1 provides the standard 
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Table Al. Correlation matrix of the AzTEC/SHADES differen- 
tial count bins as given in Table l4l The last row gives the standard 
deviation of each bin as determined through the bootstrap sam- 
pling method of Section 14,11 Together, these values can be used 
to create a covariance matrix; however, we note that the prob- 
ability distribution is not strictly Gaussian, particularly for the 
less populated bins, as evidenced by the asymmetric uncertainty 
intervals of Table [4] 
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deviation of the differential counts of each flux bin, as esti- 
mated in the bootstrap sampling method. These values can 
be applied to the correlation matrix to produce a covariance 
matrix for the data. However, the standard deviation is not 
an ideal representation of the true uncertainty distribution 
(Table Q due to the finite sampling of each bin, which re- 
sults in an asymmetric multinomial probability distribution 
(i.e. non- Gaussian). As discussed in Section [4.31 care must 
be taken when attempting to use this covariance matrix in 
typical fitting algorithms due to the high level of degeneracy 
amongst the bins. Larger flux bins could be used to reduce 
the bin-to-bin correlations; however, significantly larger bins 
would make flux resolution the limiting factor (with respect 
to the precision of the AzTEC/SHADES differential number 
counts estimate) for most practical applications of the data. 



